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Abstract. Differential Calculus is a staple of the college mathematics major's 
diet. Eventually one becomes tired of the same routine, and wishes for a more 
diverse meal. The college math major may seek to generalize applications of 
the derivative that involve functions of more than one variable, and thus enjoy 
a course on Multivariate Calculus. We serve this article as a culinary guide 
to differentiating and integrating functions of more than one variable - using 
differential forms which are the basis for de Rham Cohomology. 



Introduction 

Differential Calculus is a staple of the eoUege matliematics major's diet. It is 
relatively easy to explain the Fundamental Theorem of Calculus: Given a differen- 
tiable function / : / ^ R defined on closed interval / = [a, b], there is the identity 

j^[f{t)]dt^f{b)~f{a). 

Eventually one becomes tired of the same routine, and wishes for a more diverse 
meal. The college math major may seek to generalize applications of the deriva- 
tive that involve functions of more than one variable, and thus enjoy a course on 
Multivariate Calculus. Actually, there is a veritable buffet of ways to differentiate 
and integrate a function of more than one variable: there is the gradient, curl, 
divergence, path integrals, surface integrals, and volume integrals. Plus, there are 
many "Fundamental" Theorems of Multivariate Calculus, such as Stokes' Theorem, 
Green's Theorem, and Gauss' Theorem. 

We serve this article as a culinary guide to differentiating and integrating func- 
tions of more than one variable - using differential forms which are the basis for de 
Rham Cohomology. 

Gradient, Curl, and Divergence 

Let's focus on functions of three variables. First, let's fix a "simply connected" 
closed subset Z? C K^; this means we can define an integral between two points 
without concern of a choice of path in D which connects them. (For example, a 
subset in the form D = [a, b] x [c, d] x [p, q] is simply connected.) Recall that the 
gradient of a scalar function / : _D ^ R is the vector- valued function V/ : D ^ M.^ 



defined by 



ox ay oz 



(Here i = (1, 0, 0), j = (0, 1, 0), and k = (0, 0, 1) are the standard basis vectors for 
3-dimensional space.) The curl of a vector field F : £) — > R'^, which we write in the 
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form F 

by 



Mi + A^j + Pk, is the vector- valued function V x F : £) 



defined 
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(If you don't remember how to compute 3x3 determinants, don't worry; we have 
a short exposition on them contained in the appendix.) Finally, the divergence of 
a vector field G : Z? ^ M"^, which we write in the form G = 5i + Tj + ?7k, is the 
scalar function V ■ G : D ^ R defined by 

V G- — + — + — 

dx dy dz 

Larson j4] provides more information on Multivariable Calculus. Marsden and 
Tromba [5] is also an excellent source. Hughes-Hallet et al. [3] provides a novel 
approach to the subject using concept-based learning. 



We would like to answer the following questions: 

(1) Which F : ^ have a scalar potential f : D such that F = Vf? 

(2) Which G : D ^ have a vector potential F : D ^ M.^ such that G = 
V x F? 

(3) Which f : D ^R have a vector potential G : D ^R^ such that / = V ■ G ? 

Relating Gradients and Curls. Let's answer the first of our motivating ques- 
tions. We'll show the following: 



F — \7 f is a gradient if and only if the curl V x F = 0. 

First assume that F = V/. That is, F = M i iVj Pk = ^ i + ^ j + ^ k. 

dx dy dz 

To find the curl of F, we have 

VxF-P^ 9iV\. /dM dP\. f dN dM 
dy dz ) \ dz dx ) ^ \dx dy 
d^f d^f \. f d^f d^f \ . f d^f d'^f 



dy dz dz dy J \dz dx dx dz J \ dx dy dy dx 
= 0. 

Note that we can change the order of differentiation as long as / is continu- 
ously twice-differentiable: with this assumption, Clairaut's Theorem states that 
the mixed partial derivatives are equal. (This theorem and its proof can be found 
in pgs. 885 and A-48].) 

Now what if V X F = 0? Can we find an / such that F = V/? Well, if V x F = 0, 
then 

'dP dN\_/dM dP^^ fdN dM\ _ ^ 
dy dz J \ dz dx J \ dx dy 
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Now define the scalar function f = u + v + 'w\'ii terms of the definite integrals 

M{a, y, z) da, 



u{x,y,z) 
v{x,y,z) 
w{x,y,z) 



N(X,T,Z} - —(X,T,Z) 

dy 



du dv 

P{x,y,Q - -7r{x,yX) - -^{x,yX) 
oz dz 



dC 



where we have fixed a point a;o i + J/o j + •^o k G D. (Here's where we subtly use the 
assumption that D is simply connected: the integrals are independent of a choice 
of path in D which connects xq i + yo j + 2:0 k and xi + y] + zh.) Upon changing 
the order of differentiation and integration, we have the derivatives 
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(We can only do this if M, N, and P 


are 



dP dN 
dy dz 



also follows from Clairaut's Theorem.) It follows that F ~ V/: 



dC = 0. 











du dv dw 
dz dz dz 



du dp dw . Idu dv d 
dx px /OX J \_dy dy /Jy 

We conclude that F : D — > M'^ has a scalar potential / : D 
F V/ if and only if the curl V x F = 0. 



k = Mi + iVj + Pk. 

R such that 



Example. Let's see this result in action. Consider the scalar function 



f{x,y,z) = y/x^ +y^ + z2. 
We compute the gradient using the Chain Rule: 



9/ 

dx 

df 



{x^ +y^ + z^) 



-1/2 



2 



2x = 



dz 2 
This gives the vector field 



dy 

df 1/21 2 I 2\ 
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• 2z = 



X 




+ y'^ + z'^ 
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F = V/ = 



xi + yj 



x'^ + y'^ + z"^ 

(See Figure [T] for a plot. To keep with the gastronomical motif of this article, 
perhaps we should call this direction field a "Prickly Pear" ?) Recall that the curl 
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Figure 1. Plot of F = V/ for /(x, y, z) ^ + y'^ + z'^ 
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V X F 
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Note that the other mixed partial derivatives give rise to a similar function. Thus, 
V X F 0. 



Relating Curls and Divergence. Let's return to the second of our motivating 
questions. We'll show the following: 
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G = V X F is a curl if and only if the divergence V • G = 0. 



First assume that G = V x F for some F = Mi + Nj + Pk. That is, 





fdM 






[ dz 


^ dx 



dx dy ) 



We compute the divergence as 

_ dS dT dU 
ox ay oz 

\ dx dy dx dz 
= 0. 
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dy dz dy dx 



d^N d^M 



dz dx dz dy 



(Here, we assume that M, N, and P are continuously twice-differentiable so that 
we can interchange the order of differentiation.) 



Now what if V • G = 0? In this case, we write 
dS dT dU 



dx dy dz 



= 



Uix,y,z) = U{x,y,zo) 





'dS 


dT' 


f 


dx 


dy 



for any fixed point xq'i + yoj + zoii G D. Define the function F — M i + N j + Ph 
in terms of the definite integrals 

rv 



M{x, y,z)^ / T{x, y, C) dC, 



U{x, T, zq) dr 



Vo 

N{x,y,z) = - / S{x,y,C)dC 

J Zo 

P{x,y,z) = 0. 

Upon changing the order of differentiation and integration, we have the derivatives 
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It follows that V X F = G. 



dP 


dN 


dy 


dz 


dM 


dP 


dz 


dx 


dN 


dM 


dx 


dy 





'dS 


dT' 


f 

1 zo 


dx 


dy _ 



dC = U{x,y,z). 



Exercise. Fix three real numbers a, 6, and c such that a + 6 + c = 0, and define 
the function 

G = aa;i + 6yj + czk. 

(See Figure [2] for a plot.) Show that V • G = 0. Moreover, find a function F such 
that G = V X F. {Hint: Guess a function in the form F = Ayzi + Exzi + CxyVi 
for some real numbers A, B, and C.) How many such functions F do you think 
there are? 
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Figure 2. Plot of G = xi + 22/j - 3zk 




Example. Consider the vector function 
a; i + y j + z k 



in terms of f{x, y, z) — \/ + y"^ + z'^ 



-\- y"^ -\- 

(Remember the "Prickly Pear" ?) We will show that there is no vector potential F 
such that G = V X F. The idea is to show that the divergence of G is nonzero. To 
this end, we compute the higher-order partial derivatives using the Quotient Rule: 
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(Check as an exercise!) Hence we find the divergence 

V C - I I - 2^' + 2?/' + 2z" _ 2 

9x2 ay2 a2;2 (3.2 _^ y 2 _^ ,2)3/2 J' 

Since V • G 7^ 0, there cannot exist a function F such that G = V x F. 
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Differential Forms 

Naturally, the three motivating questions in the introduction are rather naive 
questions to ask ~ although a bit difficult to answer! - so let's make them a little 
more interesting. We'll rephrase the definitions above in a way using differential 
fc-forms; this will make integration more natural. 

• A 0-form is a continuously differentiable function f : D ^ M.. These are 
the functions we know and love from the one- variable Differential Calculus. 

• A 1-form is an arc length differential in the form rj = M dx + N dy + P dz 
for some continuously differentiable vector field F : Z3 — > R'^ in the form 
F = Mi + A^j + Pk. Note that we can express this in the form rj — F ■ dr 
using the dot product. We'll use this to define path integrals later. 

• A 2- form is an area differential in the form uj — S dy dz -\-T dx dz + U dx dy 
for some continuously differentiable vector field G : D ^ R'^ in the form 
G = S i + T j + U k. Note that we can express this in the form cu = G ■ dA. 
As the notation suggests, we'll use this to define surface integrals. 

• 3-form is a volume differential in the form v = f dx dy dz for some contin- 
uously differentiable function / : R. Note that we can express this 
in the form ly — fdV. As the notation suggests, we'll use this to define 
volume integrals. 

We'll denote n''(D) as the collection of fc-forms on D for fc = 0, 1, 2, and 3. 
Boothby [2] provides a rigorous treatment of differential forms. 

We claim that each il''{D) these is a linear vector space. This means that the 
linear combination of two /c-forms is another /c-form. Consider, for example, 1-forms 
and 2-forms. Given scalars a and f3, we have the identities 

a {Ml dx + Nidy + Pi dz) + (3 {M2 dx + N2dy + P2 dz) 

= {a Ml + (3 M2) dx+{aNi+p N2) dy + {a Pi + (3 P2) dy 

a {Si dy dz + Ti dx dz + Ui dx dy) + (3 {S2 dy dz + T2 dx dz + U2 dx dy) 

= {a Si + /3 S2) dy dz + {aTi + (3 T2) dx dz + {a Ui + /3 U2) dx dy 

Thus ^l''{D) is a linear vector space for k = 1,2. The argument for fc = 0,3 is 
similar. 

Maps Between DifTerential Forms. How are these linear vector spaces related? 
Well, there are linear maps between them! Consider the following "differential" map 
d : n°{D) n^{D) from 0-forms to 1-forms: 



/ 

I 

Of , df 
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This is hnear because d(a fi + fi f2) — a dfi + (5df2- Here, we write dr — dxi.+ 
dy] + dzk as the arc length differential - expressed as a vector. Hence we find the 
answer to our first motivating question: 

The vector-valued function F = V/ has a scalar potential f if and only if the 
l-form F ■ dr = df is the differential of some 0-form f. 

Similarly, consider the following "differential" map d : n^{D) n^{D) from 
1-forms to 2-forms: 

7] = Mdx + Ndy + Pdz = F -dr 

i 

fdP dN\ , , fdM dP\ , , fdN dM\ , , 
= (V X F) • dA. 

This is linear because (i(a ?7i +/3 772) = adr}\+ fi dri2- Here, we write dA = dydzi + 
dxdz} + dx dy k as the area differential - also expressed as a vector. Hence we find 
the answer to our second motivating question: 

The vector-valued function G = V x F has a vector potential F if and only if the 

2- form G • dA = dr] is the differential of some l-form rj = F ■ dr. 

Recall that F = V/ is a gradient if and only if the curl V x F = 0. In other words, 
given a l-form r] = M dx -\- N dy -\- P dz = F ■ dr, we have r/ = df as the differential 
of a 0-form if and only if the 2-form drj = 0. This is yet another answer to our first 
question! 

Finally, consider the following "differential" map d : fl'^{D) — > fl^{D) from 2- 
forms to 3-forms: 

LO = S dydz -\-T dxdz -\-U dxdy = G ■ dA 

i 

(dS dT dU\ , , , 
duj=\— + — + —] dxdydz = ( V ■ G) dV. 

\ox oy oz J 

Again, this is linear because d{auji + (3uj2) = adioi + (3djjj2- Remember that 
dV = dxdydz is the volume differential. Hence we find the answer to the last of 
our three motivating questions: 

The vector-valued function / = V ■ G has a vector potential G if and only if the 

3- form f dV = du is the differential of some 2-form co = G ■ dA. 

Recall that G = V x F is a curl if and only if the divergence V • G = 0. In other 
words, given a 2-form cj = S dy dz -\- T dx dz -\- U dx dy = G ■ dA, we have w = dri 
as the differential of a l-form if and only if the 3-form doj = 0. This is yet another 
answer to our second question! 
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We'll summarize this via the following diagram of maps: 

{0} . n^{D) n'{D) n^D) ^^isie^ ^^^d) > {o}. 

All of the various definitions of differentiation for functions in three variables are 
all related via linear transformations on the vector spaces of fc-forms! 

Example. Let's consider the differential 1-form 

ri = yzdx + xzdy + xy dz. 

This is in the form M dx + N dy + P dz for M = y z, N = x z, and P = xy. (See 
Figure [3] for a plot of F = Mi + iVj + Pk.) We'll compute the differential dry. We 
have the partial derivatives 

dP ON dM dP ON dM 

"a H-=2;-x = 0, — — =y-y = 0, — — =z-z = 0. 

oy oz az ox ox ay 

Hence dij = 0. 



Figure 3. Plot ofF = y zi + x zj + xyk 




Exercise. Fix three real numbers a, 6, and c such that a+b+c = 0, and differential 
2-form 

ll> = ax dy dz + by dx dz + c z dx dy. 
Show that Lu = drj is the differential of the 1-form ri = byzdx — axzdy. Is this the 
only 1-form rj such that uj — drj? 
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Integration 

We know now how the various forms of differentiation are aU related, but what 
about integration? First, we review what we mean by integration. As I? C M'' is 
a subset of three-dimensional space, we can define integration for one variable, two 
variables, or even three variables. 



Fundamental Theorem of Calculus. In order to define an integral for one vari- 
able, let r : / ^ Z) be a continuously differentiable map, which we write in the form 
r — xi + j/j + zk, defined on a closed interval / = [a,b]. We say that the image 
C C D is a path. Given vector field F : D ^ M"^, the composition gives the 1-form 
r] = F(r(t)) • dr, so naturally a path integral as 



c 



Mm I 



dt. 



Recall that F V/ has a scalar potential / if and only if jj 
of a 0-form /. In this case, the path integral simplifies to 



df is the differential 



d_ 

dt 



f{r{t)) 



dt 



df = /(r) 



c 



dC 



/(r(6))-/(r(a)). 



Hence the integral is independent of the path C, and it is only dependent on the 
endpoints dC — {r(a), r{b)}. This is just the Fundamental Theorem of Calculus. 



- r^ and z = 0. We may think of this as being the image of the map 
which sends t r costi -I- r sinij. The circle is just an example 



Example. Let _D = be all of three-dimensional space, and let C denote the 
circle of radius r in the plane. Note that C is the set of points xi + yj + zk such 
that -I- 
r : [0, 27r] ^ 
of a path! 

Now let's consider a vector field F = — (j//2) i+(a::/2) j. (See Figure|4]for a plot.) 
We'll explain why there does not exist a scalar function / such that F = V/. The 
idea is to suppose that / does indeed exist and then compute the integral of F 
around the path C. Let's consider the following differential 1-form: 



r] = F(r(i)) • dr 



r cos t . r sin t , . . 
— 1 -I — J ] ■ [ —r cos tdti + r sm < dt j 



cos^ tdt + r"^ sin^ t dt r 



= — dt. 
2 



If F = V/, then r; = /(r(27r)) - /(r(0)) = 0. But actually. 



2^ ^2 

c Jo ^ 



dt = TT r"^ . 



Hence there is no function / such that F = V/. (Of course, we could have seen 
this sooner by computing the curl V x F = k and realizing it's a nonzero vector.) 



Stokes' Theorem and Green's Theorem. In order to define an integral for two 
variables, let r : i? ^ _D be a continuously differentiable map, which we write in 
the form r = a;i-|-yj-|-2:k, that is defined on an closed region R = [a, h] x [c, d] . We 
say that the image S' C D is a surface, and assume that its boundary C = dS is a. 
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Figure 4. Plot of F = -(y/2) i + (a;/2) j 




du dv. 



curve as before. Given vector field G : Z) ^ M"^, the composition gives the 2-form 
u! — G(r{u,v)^ ■ dA, so naturally surface integral is defined as 

S{r{u,v)) r(r(M,u)) U{r{u,v)) 

dx dy dz 

du du du 

dx dy dz 

dv dv dv 

Recall that G = V x F has a vector potential F if and only ii u; = drj is the 
differential of a 1-form rj. In this case, the surface integral simplifies to 

(V X F) • dA = lldr]= 1 V= [ F ■ dr. 

IS JJs JdS Jc 

Hence the integral is independent of the surface S, and it is only dependent on the 
boundary C — dS. This is known as Stokes ' Theorem. 

Let's consider a special case, where r : i? ^ actually maps into the plane. 
Then S* C K^, so write Y — M\ + Nj. Stokes' Theorem in the plane reduces to the 
statement 

'dN dM\ 



, dxdy = di]^ ?]= [Mdx 
s \dx dy J JJs JdR J as 



Ndy]. 



This is known as Green's Theorem. Similarly, for the orthogonal vector field F^ = 
— iVi + M j, the Divergence Theorem in the Plane is the expression 

Example. Let = R"^ be all of three-dimensional space, and let S denote the 
disk of radius r in the plane. The latter is just the set of points x i + y j + z k such 
that x^ + y^ < and z — Q. We may think of this as being the image of the map 
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9 

y =r 



p coaOi + p sin 6]. Note that the 
^ which wc considered before. 



r : [0, r] x [0, 27r] which sends (p, 9) 

boundary C = dS is simply the circle x'^ 

Let's return to the vector field F = — (y/2) i + {x/2)j. We'll consider the dif- 
ferential 2-form w = (V X F) • dA = k • dA = dx dy. Stokes' Theorem (or really 
Green's Theorem, since it's in the plane) states that 



lis Ic 



7rr 



F-dr 

IS Jc 

(Recall the previous example.) Let's compute the integral on the left-hand side in a 
different way. This differential 2-form can also expressed in polar coordinates using 
the determinant 

T{rip,0)) U{v{p,9)) 

dy dz 
dp dp 
dy dz 

~ de 



s{v{p,e)) 

dx 

d~p 

dx 

de 



de 



dpde 





cos^ 
-p sin e 





sin^ 
p cose 



dpde 



= pdp de. 

Hence wc have the integral 



27r 



pdpde 



pdp 



271- 



de 



TTr 



dx dy = 
IS JJS Jo Jo 

Of course, this is just the area of the disk S. 

Exercise. Let S be any surface in the plane with a boundary C = dS. Show 
that its area can be computed using a path integral around along C. That is, 

/■ —ydx + xdy 

Ic 2 ■ 



Area(S') 



dx dy 



Gauss' Theorem and the Divergence Theorem. In order to define an integral 
for three variables, let r : B — > be a continuously differentiable map, which we 
write in the form r = xi + yj + zk, defined on an closed region B = [a,b] x 
[c, c/] X \p,q\. We say that the image D C is a region, and assume that its 
boundary S = dD is a surface as before. Given a scalar function / : — > M, the 
composition gives the differential 3-form i/ = /(r(u, w,w)) dV, so naturally define 
a volume integral as 



D 



f{r{u,v,w)) 



dx 


dy_ 


dz 


du 


du 


du 


dx 


dy 


dz 


dv 


dv 


dv 


dx 


dy 


dz 


dw 


dw 


dw 



dudvdw. 



(As with surface integrals, we've expressed the integrand above using a 3 x 3 deter- 
minant. Most texts refer to this as the Jacobian of the transformation r : B — > D.) 
Recall that / = V • G has a vector potential G if and only if z/ = cL; is a differential 
of a 2-form u). In this case, the volume integral simplifies to 

/// {V-G)dV= [[[ du)= [[ w= [[ G-dA. 
JJJd JJJd JJdD JJs 
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Hence the integral is independent of the region V, and it is only dependent on the 
boundary S = dD. This is known as Gauss' Theorem, or the Divergence Theorem. 

Example. Let D denote the solid sphere of radius r in three-dimensional space. 



This is just the set of points x i + y j + z k such that x'^ + + z'^ < . 
think of this as being the image of the map r : [0, r] x [0, tt] x [0, 2it] 
sends 

{p, (l),6) ^ p sin (j) cos 6i-\- p svacj) sm9] + p cos 4> k. 

Note that the boundary S = dD is simply the sphere x"^ + + z"^ = r^. 
Let's return to the vector field 
xi + yj + zk 



We may 
R3 which 



G = Vg = 



in terms of g{x, y, z) = \/x'^ + y"^ + . 



\fx?^^r]p^A-Z^ 

(The "Prickly Pear" again!) Wo will compute integrals and verify the Divergence 
Theorem. We may express the volume element dV = dx dy dz in spherical coordi- 
nates using the determinant 



dV 
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de 



dp d(p de 



coau sm(p 
p cos e cos 4 
—p sine sine 



smy smcp 
p sin e cos ( 
p cos e sin < 



-p smc 




dp d<p de 



= p sin (j) dp dip d9. 
Hence we have the differential 3-form 

v = \7 ■ G{r{p,(l),e)) dy {P^ sin <j) dp d<j)de) 

(Recall that V ■ G = 2/g.) This gives the integral 



2 p sin (j) dp d(p d9 . 
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Jo 



2 p sin (/> dp d(l) dO 











- /■27r 


I 2pdp 




/ sin (j) dtj) 




/ de 


Jo 




Jo 




Jo 



= 4 TT r . 

On the other hand, we may express the area element is the determinant 



dA 



i 


j 


k 


dx 


dy 


dz 


d4 
dx 


d(f> 
dy 


d4> 
dz 


'de 


de 


de 



d(j)de 



p cos e cos (f. 
—p sin 9 sin c 



p sm e cos ( 
p cos e sin ( 



-p sine 




— p sin cf)[p sin (p cos ei + p sin (j) sin j + p cos ((> k) 



^de. 



Hence we have the differential 2-form 

uj = G (r{p, 4>,9)) ■ dA = - ■ [p sin(j)r de ( 

This gives the integral 

// w = I I sin(l)d9d(j) = 

J JdD Jo Jo 



p^ sin 



d9. 



d(l)de 







- ^27r 




/ sin (j) d(j) 




/ de 


= 47rr^ 


Jo 




Jo 
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(Remember that p = r along the boundary S — dD.) This indeed verifies that 

N-G)dV= ff G-dA. 



D 

Arfken and Weber [1] provides a plethora of formulas for integrating in different 
coordinate systems other than spherical. 

Exercise. Let D be any region in M.^ with a boundary S — dD. Show that its 
volume can be computed using a surface integral on S. That is, 



Vol(L») = / / / dxdydz 



D 



X dy dz -\- y dx dz + z dx dy 



The Moral of the Story 

In this paper we have shown that all of the theorems - those for differentiation 
and integration - can be expressed using differential forms. Here's the idea in a 
nutshell for any region D. We consider a series of linear vector spaces fl''{D), where 
we have "differential" maps 

... , n''-^{D) n'^iD) n'^+^D) > ••• 

We would like to know the answer to the following question involving differentiation: 

Which k-forms uj G ^'^'^^{D) are in the form 
iu = drj for some {k — l)-form rj £ il''{D) ? 

The partial answer should be: 

If UJ = drj for some (k — l)-form rj, 
then du> = Q as a {k + \)-form. 

A complete answer involves computing something called de Rham Cohomology. 
(The desert after the Main Course?) For instance, we have computed de Rham 
Cohomlogy in this article for certain subsets of three-dimensional space. We expect 
to generalize the integration formulas above by saying something like 

/ w = / d77 = / 77 

Jd JD JdD 

SO that the integral would be independent of the region D, and would be only 
dependent on the boundary dD. This is the Generalized Stokes' Theorem. Surpris- 
ingly, this entire theory can be worked out for "many" sets D. But don't take our 
words for it - just a friendly Differential Geometer! 

Appendix: 3x3 Determinants 
A 3 X 3 determinant is the quantity 

aii(a22 - a32) - ai2(a23 031 - 021 033) 
+ ai3(a2i a32 - a22 asi)- 



Oil ai2 ai3 
021 a22 a23 
031 a32 a33 
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Another way is to compute the determinant is by using the following diagram: 




an ai2 ais an ai2 




XX/ 



a21 a22 023 «21 022 

X X X \ 

C(31 0'32 Oaa 0,31 a32 




Here we multiply along the arrows, then add or subtract depending on the direction 
of the arrow: 

(an ^22 ass + 0,12 a2s asi + ais 021 032) 

- (an 0123 0132 + ai2 021 033 + ais 022 asi) . 
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